First application of plutonium in soil erosion research on terraces 
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Abstract: The spatial distributions of 73°**4°Pu and '°’Cs in soils from Longji Rice Terraces were investigated to evaluate soil 
erosion. The activity concentrations of 7°'*4°Pu and '°’Cs in the surface soils of the paddy fields were in the range of 0.089-0.734 
and 1.80-7.88 mBq/g, respectively. The activities of 7°'**°Pu and '°’Cs showed very similar distribution trends, first increasing and 
then decreasing with increasing elevation. The **°Pu/**Pu atom ratios in the surface soils ranged from 0.162—0.232. The activities 
of 73°*74°Py and !°’Cs in the soil cores tended to be uniformly distributed within the plowed layer and declined exponentially below 
this depth. The mean soil erosion rates of Longji Rice Terraces estimated by 7°°*°Pu and '°’Cs tracer methods were 5.44 t/(ha-a) 
and 5.16 t/(ha-a), respectively, which demonstrated that plutonium can replace !3’Cs as an ideal tracer for soil erosion research in 
the future. Landform features are the main factors affecting the distribution of plutonium and !3’Cs as well as soil erosion in the 
Longji Rice Terraces. 
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1. Introduction 


Soil erosion leads to a reduction in soil fertility and degradation of cultivated land, which severely affects agricultural land. 
Globally, approximately 1.9 billion hm? of agricultural land is affected by land degradation and is increasing at a rate of 
approximately 5—7 million hm? per year [1]. Human activities and global warming have resulted in increased soil erosion. Therefore, 
it is necessary to conduct soil erosion research to provide basic data to prevent land degradation. 

Radionuclides have been proven to be powerful tracers for soil erosion investigation [2-4]. Plutonium and !37Cs (712 = 30.17 a) 
are artificial radionuclides mainly produced by artificial nuclear activities [5-9]. 7°Pu (71/2 = 2.411 x 104 a) and “Pu (712 = 6563 
a) were the most important plutonium isotopes. The *°?**°Pu released to the environment by atmospheric nuclear tests was 
approximately 1.26 x 10!° Bq [10]. The atom ratio of 7°Pu/??°Pu has become an important factor in investigating the source and 
migration processes of plutonium in the environment [11]. Accelerator mass spectrometry (AMS) is an effective and robust 
technique for studying long-lived actinides (e.g., 7°°U, 7>7Np, 7°°Pu, and *“°Pu) [12-14], and inductively coupled plasma mass 
spectrometry (ICP-MS) is currently the most widely used technique for measuring actinides because of its lower cost compared to 
AMS. The distribution, sources, migration, and applications of !3’Cs have attracted extensive attention over the past few decades 
[15-20]. 7°°'*4°Pu and !7Cs in the environment are rapidly and strongly adsorbed by fine soil particles and redistributed with the 
soil itself [21-23]. Thus, soil erosion can be traced using the inventories and distribution characteristics of 77°‘*4°Pu and !3’Cs in the 
soil [21, 24-29]. 

Although !?7Cs has been extensively applied to soil erosion in past decades [1, 24, 25, 30-32], 70% of the total !37Cs deposited 
by atmospheric nuclear tests in the environment has decayed away so far, which has significantly reduced the sensitivity of '°’Cs 
as a tracer of soil and sediment transport in the future. 7°°'*°Pu is an ideal replacement for '3’Cs due to its longer lifetime and the 
same origin and strong association with soil particles as '3’Cs. In recent years, plutonium has been successfully applied to soil 
erosion studies in the hilly landscape north of Munich, Germany [33, 34]; the Alpine grasslands of Switzerland [27]; the lower 
Cotter catchment of the Australian Capital Territory [35]; the hillslopes of the Zarivar Lake watershed, Kurdestan Province, Iran 
[24]; the coastal area of Liaodong Bay in Northeast China [4, 21, 28]; and the hilly gully region of the Loess Plateau in China [36]. 
However, the study areas are mostly limited to grasslands, woodlands, hillslopes, and commonly sloping cultivated land. To the 
best of our knowledge, the application of plutonium to soil erosion on terraces of stepped-sloping cultivated lands has not yet been 
published. Terraces are fields reclaimed on the slopes of mountainous areas where precipitation is abundant. The steepness of the 
terraces depends on the slope gradient of the original hillside. Excess water from the upper terraces flows to the lower terraces due 
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to gravity. Terraced soil is affected by tilling, loose, and easily washed and carried by rain, particularly during periods of heavy 
rainfall. Therefore, the steeply sloping rugged mountain topography, loose soil, abundant precipitation, and water flow from the top 
to bottom create favorable conditions for terrace soil erosion. Soil erosion is directly related to local crop yields and sustainable 
cultivation owing to the rugged hillside terrain, top-down diffuse water flow, and primitive coupled tillage of the Longji Rice 
Terraces; therefore, it is necessary to explore a proven means of monitoring local soil erosion. The combination of the long-lived 
artificial radionuclides 7*°*?4°Pu and '°’Cs for studying the long-term erosion of terraces would be more convincing. In this study, 
the activity concentrations of °**4°Pu and '°’Cs, and the atomic ratios of **°Pu/*3°Pu in soils from the Longji Rice Terraces in 
Guangxi, China, were measured for the first time. This study analyzed the sources and distribution characteristics of 73°*”°Pu and 
137Cs in the soils of the Longji Rice Terraces, evaluated the soil erosion status, and verify the feasibility of 73°*74°Pu instead of '°’Cs 
to study the long-term erosion of the terraces by comparing the erosion rates obtained from the °**4°Pu and !37Cs tracer methods. 


2. Material and methods 


2.1. Study area 


Longji Rice Terraces (109°32'-110°14' E, 25°35'-26°17' N) is located in Longji Town, northwest of Guilin City, Guangxi, China. 
It is a field developed on Longji Mountain, which belongs to the Yuecheng Mountain range and is the northwest branch of the 
Nanling Mountain system. The Longji Rice Terraces mainly include Longji Guzhuang Village terraced fields (LGVTF), Ping’an 
Zhuang terraced fields (PZTF), and Jinkeng Red Yao terraced fields (JRYTF). Longji Rice Terraces are in a subtropical monsoon 
climate zone affected by the convergence of cold air from the north and warm and humid air from the south, with more frontal rain 
and fog weather and an annual precipitation was 1800—1900 mm. The elevation of the Longji Rice Terraces is 300—1100 m. Longji 
Mountain is covered with lush vegetation, so it can store water well. The upper layer of Longji Mountain is permeable sandy loam, 
which can store water, and the lower layer is impermeable sedimentary rocks, shale bedrock, and slate with good water 
impermeability, located approximately 1.5—2 m below the surface soil. The excess rainwater that infiltrated the sandy loam flowed 
out of the cracks in the mountain and down to each terrace. The farming methods of the Longji Rice Terraces are primitive pair 
plowing and traditional flood irrigation. 


2.2. Sample collection 


Soil samples were collected from Longji Rice Terrace in November 2021. The sampling locations and information are shown in 
Fig. | and Table |, respectively. The sample IDs, beginning as “P” and “D”, were taken from PZTF and JRYTF, respectively. The 
sampling 


represent sites of surface soils and cores, respectively. The red lines represent the sampling sites along the slope direction. (a) PZTF. The yellow line 
represents the contour line at an elevation of 860 m. (b) JRYTF. The yellow line represents the contour line at an elevation of 970 m. 
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sites on the PZTF (open triangles and pentagons) and JRYTF (solid triangles) are shown in Fig. 1. Triangles and pentagons represent 
the surface soil and core sites, respectively. Red and yellow lines represent the sampling sites along the slope direction and contour 
lines, respectively. Fourteen paddy field surface soil samples (0—5 cm, marked by open triangles) and two paddy field soil cores 
(P1#1 with a depth of 50 cm and P1#2 with a depth of 40 cm, marked by open pentagons) were collected from the PZTF, and 19 
paddy field surface soil samples (0-5 cm, marked by solid triangles)) were collected from the JRYTF. Surface soil samples were 
collected along a specified slope direction. Samples from the PZTF and JRYTF were collected from high and low elevations along 
the southwest and southeast slopes, respectively. All soil cores were sliced at 5 cm intervals. 


Table 1 General information of the sampling sites and results of Pu and !3’Cs in soil samples. 


Sample ID Longitude(E) Latitude(N) Elevation 2394240 Dy BICs 40Py/?°Pu OM*(%)  Slope(°) 
(m) (mBq/g) (mBq/g) atom ratio 
Pl#1_0-5 110° 7’ 29" 25° 45'39" 930 0.207 + 0.006 4.04 + 0.20 0.197 + 0.015 9.3 - 
Pll 110° 7’ 29" 25° 45'37" 920 0.245 + 0.009 6.16 +0.25 0.189 + 0.015 9.1 9.4 
P12 110° 7' 28” 25° 45'36" 900 0.298 + 0.010 6.35 + 0.25 0.203 + 0.018 9.4 25.1 
P13 110° 7'27" 25° 45'33" 870 0.335 + 0.011 7.88 + 0.28 0.194 + 0.015 11.3 17.5 
PH11 110° 7’ 27" 25° 45'32" 860 0.330 + 0.012 6.20 + 0.25 0.181 + 0.015 9.9 18.4 
P14 110° 7' 26” 25° 45'28" 840 0.259 + 0.009 5.22 + 0.23 0.173 + 0.012 8.9 9.2 
P15 110° 7' 26” 25° 45'24" 810 0.228 + 0.006 4.35+0.21 0.193 + 0.011 11.1 14.0 
P16 110° 7’ 26" 25° 45'23" 790 0.089 + 0.006 1.80 £0.13 0.208 + 0.023 9.9 33.6 
PH12 110° 7' 29” 25° 45'32" 860 0.335 + 0.009 6.31 £0.25 0.192 + 0.013 11.1 - 
PH13 110° 7’ 29" 25° 45'32" 860 0.249 + 0.007 5.20 + 0.23 0.196 + 0.013 10.3 - 
PH14 110° 7’ 26" 25° 45'32" 860 0.265 + 0.009 3.40 40.18 0.187 + 0.010 10.5 - 
PH15 110° 7' 26” 25° 45'32" 860 0.264 + 0.006 6.27 + 0.25 0.196 + 0.010 13.0 - 
PH16 110° 7’ 27" 25° 45'36" 860 0.232 + 0.011 4.98 + 0.22 0.209 + 0.022 14.4 - 
P1#2_0-5 110° 7'7" 25° 45'32" 860 0.215 + 0.016 5.13 + 0.23 0.201 + 0.018 8.0 - 
D1#1_0-5 110° 9'1” 25° 48'58" 1080 0.224 + 0.012 4.99 + 0.22 0.216 + 0.026 10.0 - 
DA8 110° 9' 6” 25° 49'4" 1070 0.089 + 0.007 2.16 +0.14 0.181 + 0.029 9.0 2.4 
Dil 110° 9'2” 25° 48'56" 1050 0.105 + 0.008 2.23 £0.15 0.197 + 0.027 9.5 4.3 
DAT 110° 8' 54” 25° 4855" 1030 0.158 + 0.010 3.24 + 0.17 0.185 + 0.028 9.0 4.7 
D12 110° 9'3” 25° 48'55" 1020 0.161 + 0.012 3.79 +0.19 0.232 + 0.033 8.6 2.1 
DA6 110° 8' 54” 25° 48'51" 990 0.157 + 0.010 4.49 + 0.21 0.188 + 0.031 8.8 5.8 
D13 110° 9' 6” 25° 48'52" 980 0.162 + 0.014 4.39 + 0.21 0.204 + 0.031 8.1 1.6 
D1#2_0-5 110° 9' 8” 25° 48'51" 970 0.214 + 0.017 4.92 + 0.22 0.176 + 0.034 9.3 8.5 
D14 110° 9'7” 25° 48'48" 940 0.247 + 0.015 5.96 + 0.25 0.201 + 0.031 9.3 17.5 
DA2 110° 9' 6” 25° 48'48" 905 0.243 + 0.010 6.14 + 0.25 0.205 + 0.021 8.7 49.3 
D15 110° 9'7” 25° 48'43" 890 0.318 + 0.024 6.04 + 0.25 0.207 + 0.039 9.3 5.6 
D16 110° 9' 9” 25° 48'38" 830 0.190 + 0.011 4.94 + 0.22 0.191 + 0.024 9.2 20.3 
D1#3_0-5 110° 9' 9” 25° 48'37" 810 0.109 + 0.010 3.06 + 0.17 0.188 + 0.029 8.8 33.6 
DH11 110° 9'8” 25° 48'52" 970 0.347 + 0.023 3.44 +0.18 0.162 + 0.013 7.0 - 
DH12 110° 9'7” 25° 4851" 970 0.380 + 0.053 4.55 +0.21 0.197 + 0.025 9.8 - 
DH13 110° 9' 7" 25° 48'50" 970 0.563 + 0.057 4.29 + 0.20 0.228 + 0.032 9.9 - 
DH14 110° 9'5” 25° 48'51" 970 0.324 + 0.030 3.86 40.19 0.192 + 0.021 10.2 - 
DH15 110° 9' 4” 25° 4851" 970 0.464 + 0.028 3.87 40.19 0.189 + 0.015 9.0 - 
DH16 110° 9' 4” 25° 48'51" 970 0.734 + 0.046 4.36+0.21 0.203 + 0.022 9.0 - 


è Organic matter content. 


2.3. Sample preparation and measurement 


The samples were heated at 100 °C for 12 h until dryness, then crushed and sieved. After that, the '°’Cs activities of the 50 g fine 
soil samples were measured by an y spectrometer equipped with a high-purity germanium detector and a lead shielding chamber at 
a time of 86400 s. The detailed methodology and spectrometer are described in our previous work [37-39]. Briefly, A HPGe detector 
(GEM-C5970; ORTEC, USA) with a counting efficiency of 45% and FWHM of 1.7 keV on 1.33 MeV energy of ®Co was used. 
137Cs content was determined by detecting the 661.62 keV signature y-ray. The detection limit for °’Cs was 0.04 mBq/g. 1 g of 
pretreated soil sample was calcined at 550 °C for 12 h to determine the organic matter content. The organic matter content was 
calculated as the mass loss before and after the ashing procedure, and was expressed as the percentage loss on ashing for the soil 
sample. 

2.5 g soil sample was weighed for plutonium isotope analysis. Details of the chemical separation and purification processes have 
been described in previous studies [37-39]. In brief, 1.82 pg **’Pu was added as a chemical yield tracer to monitor and calculate the 
chemical recovery of the plutonium isotope analysis, the main analysis processes including digestion, filtration, adjustment of 
acidity and valence, purification of plutonium with AG1-X8 anion exchange resin and re-purification of plutonium with AG MP- 
1M anion exchange resin. Finally, the concentrated plutonium solution was placed in a 1.5mL sample bottle with 4% nitric acid and 
carried to the College of Radiology and Protection, Soochow University, for measuring °°*”4°Pu activity and *“°Pu/?*°Pu atom ratio 
using sector-field ICP-MS (Finnigan Element 2, Bremen, Germany). The chemical recovery of Pu ranged from 60% to 80%, and 
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the limits of detection for ICP-MS calculated from 3 times the standard deviation of the blank solution (4% HNO3) were 0.09 fg/mL 
for >°Pu and 0.14 fg/mL for °Pu with a counting time of 150 s. The concentration of plutonium isotopes was calculated by 
subtracting the blank from the measured results of samples, and the measured signal intensities of ™°Pu and ™4°Pu in procedure 
blank were 8.5 cps and 6.4 cps, respectively. 

3. Results and discussion 


3.1. Activities of *°'“°Pu and "Cs in soil samples 


3.1.1. 7Pu/*°Pu atom ratios, activities of 7°°*°Pu and !3’Cs in surface soil samples 


The mean atom ratio of *°Pu/*??Pu in the surface soil samples was 0.195+0.014 (n=33), ranging from 0.162+0.013 to 
0.232+0.033 (Table 1), which agreed well with the value of the global fallout (0.178+0.023 (0-30°N)) [40]. This result is consistent 
with those of previous reports [41-46]. The result is higher than the **°Pu/*?°Pu atom ratio for nuclear weapons-grade plutonium 
(0.01-0.07, e.g., 0.028-0.037 for the Nagasaki atomic bomb explosion) [47, 48], and lower than those for the Chernobyl nuclear 
accident (0.403-0.412) [49] and Fukushima nuclear accident (0.303-0.330) [9], indicating there is no regional radioactive 
contamination in the studied area, so all the Pu and Cs are sourced from the global fallout, then they can be used as tracers for soil 
erosion research. 

The activities of 7°°74Pu and '°’Cs in surface soil samples are listed in Table 1. It can be seen that the activities of 7°'?*°Pu 
range from 0.089+40.007 to 0.734+0.046 mBq/g, with a mean value of 0.264 + 0.132 mBq/g (n=33), similar to that of surface soils 
(0.365 + 0.014 - 0.819 + 0.066 mBq/g) from Shaodong of Hunan Province in China [43], and within the data range (0.005-1.988 
mBq/g) of Chinese surface (< 10 cm) soil [50]. 77°*4°Pu activities in the surface soil of PZTF (0.089 + 0.006-0.335+0.011 mBq/g, 
with a mean of 0.253+0.064 mBq/g)and JRYTF(0.089+0.007-0.734+0.046 mBq/g, with a mean of 0.273+0.168 mBq/g) are similar. 
The activities of '3’Cs in the surface soil samples range from 1.80+0.13 to 7.88+0.28 mBq/g, with a mean value of 4.67+1.38 mBq/g 
(n=33) which was slightly higher than that of 3.39 mBq/g from Puding County of Guizhou Province [51]. Because all samples were 
collected from a limited spatial range with a diameter of 1000 m, the large variation in '°’Cs activity indicates local alteration after 
deposition. 
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Fig. 2 Distributions of 7°°*?4°Pu and '37Cs activity in the surface soils of paddy field from Longji Rice Terraces. (a) PZTF. (b) JRYTF. 


As shown in Fig. 2a, the activity of 73°**4°Pu in the surface soils of the PZTF first increased and then decreased from an elevation 
of 930 m to 790 m. The maximum value of PZTF appeared at elevations of 860 and 870 m. Fig. 2b indicates that 7°°**4°Pu activities 
in the surface soils of JRYTF, except for the samples at the highest elevation of 1080 m, also showed the same trend as PZTF. The 
maximum value of the JRYTF appeared at an elevation of 890 m. The same trend in !3’Cs activity is observed in Fig. 2. Similarly, 
the '°’Cs inventories of sloping farmland in the southeast of the Qinghai-Tibet Plateau showed a fluctuating increasing and 
decreasing trend along elevation [25]. All possible factors affecting plutonium and '°’Cs concentrations included soil texture, 
organic matter content, soil moisture, particle size, elevation, and slope gradient. All samples were collected from a small area of 
rice fields under the same basic environmental conditions, which had the same soil texture, soil moisture, and grain size; therefore, 
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these factors could be excluded. Elevation was also not a major factor affecting the plutonium and !37Cs concentrations in the studied 
area because the elevation span of the studied sites was within 300 m, and the concentrations of plutonium and !3’Cs did not increase 
or decrease monotonically with elevation. A strong correlation between *°*?#°Pu activity concentration and organic matter content 
was reported in our previous study [39] and by other researchers [52]. However, the surface soil samples of Longji Rice Terraces 
were all collected from rice paddies with little variation in organic matter content (Table |), suggesting that landform is more likely 
to be the main factor in 7*°*“°Pu activity variation relative to organic matter. Landform factors, such as slope gradients, are among 
the most important factors affecting the spatial variation in soil erosion rates [53, 54]. Fig. 3a shows the distributions of *°'**°Pu 
and "Cs in surface soils at an elevation of 860 m in the PZTF. The activities of °'7°Pu and !37Cs were 0.21540.016-0.33540.009 
mBq/g and 3.40+0.18-6.31+0.25 mBq/g respectively, with the mean values of 0.270+0.046 mBq/g (n=7) and 5.27+1.09 mBq/g 
(n=7). Fig. 3b shows the distributions of 7°?'74°Pu and !3’Cs in the surface soils at an elevation of 970 m in the JRYTF. The activities 
of 739*40Py and '37Cs were 0.214+0.017-0.734+0.046 mBq/g and 3.44+0.18-4.92+0.22 mBq/g respectively, with the mean values 
of 0.432+0.172 mBq/g and 4.18+0.49 mBq/g. At the same elevation and with little variation in organic matter content, the activities 
of 73°*249Py and '3’Cs at different sites showed discrepancies, indicating that the slope gradient was a significant factor in the 


discrepancy. 


0.36 ee 72394240py 0.8 ge 239+240py 
eo 137Cs 


0.34 


oad 
a 


> 

pe 

Ss 
bad 
a 


tn 


is 
` 


E 
137Cs activity(mBq/g) 


137Cs activity(mBq/g) 
2394240 py activity(mBq/g) 


2394240 bu activity(mBq/g) 
e 
: 


= 

i 

N 
ha 


0.2 


> 
i 
S 


PHI3 PHI2 PHIT = P1#2_0-5 PH14 PHIS PHI6 DHII D1#2_0-5 DH12 DH13 DÐDHI4 DHI5S ~——DHI16 
Sample ID Sample 1D 


Fig. 3 Distributions of °°*?°Pu and !?7Cs activities in the surface soils of equal elevation paddy field in Longji Rice Terraces. (a) With an elevation of 860 
m in PZTF. (b) With an elevation of 970 m in JRYTF. 


Similar to the staircase gradient, here, we define the angle (a) between the slope line of adjacent sampling sites and the horizontal 
plane as the slope gradient of the terraces between adjacent sampling sites (Fig. 4). The slope gradient was determined by the width 
(b) of the horizontal terraces at each level and the height (A) between them and the upper or lower terraces, which reflected the 
steepness of the terraced section. Fig. 5a shows the slope gradient distribution between two adjacent sampling sites in the PZTF. 
The slope gradient first decreased and then increased with elevation from 900 to 790 m, which is the opposite of the distribution 
trends shown for the activities of 73°*74°Pu and !3’Cs, indicating that the larger the slope gradient, the lower the activities of 7°°**4°Pu 
and !3’Cs and the greater the soil erosion. This is consistent with the results for sloping cultivated fields in Nanchong of Sichuan 
and Kaixian of Chongqing [55]. Khodadadi et al. [24] and Su et al. [25] reported that soil erosion affected the distribution of 
plutonium and !3’Cs in hillslope soils. Li et al. [53] and Van Oost et al. [54] also reported that slope gradient is the most important 
factor affecting the spatial variability of soil erosion rates. 


Fig. 4 Schematic diagram of the slope gradient between adjacent sampling sites of the terraces. 
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The activities of 7°°'*4°Pu and '3’Cs at an elevation of 930 m (P1#1_0-5) were lowest, except for that with an elevation of 790 m 
(P16). However, the slope gradient (9.4°) between this site P1#1_0-5 and the sampling site P11 was low (Fig. 2a and Fig. 5a). This 
is inconsistent with the rule that the greater the slope gradient, the lower are the activities of 73°*”4°Pu and !3’Cs. This may be due to 
the fact that the slope gradient between this site (P1#1_0-5, located at the top of the terrace slope closest to and directly below the 
observation platform) and the observation platform should be the largest. Serious soil erosion usually occurs at the top of slopes in 
sloping and terraced farmlands without soil material supplement [25]. This slope gradient exacerbated the scouring of the soil at 
this site by the rapid runoff above it (especially the concentrated rainwater from the observation platform), and there was insufficient 
soil replenishment at this site, resulting in the lowest activities of °°**4°Pu and '°’Cs. The slope gradient (9.2°) between the sample 
sites with elevations of 860 m (PH11) and 840 m (P14) was smaller than that (18.4°) between the samples with elevations of 860 
m and 870 m (P13), whereas the maximum activity of °*74°Pu and '!3’Cs was not at an elevation of 840 m but at the elevations of 
860 m and 870 m place. This is due to the influence of existing trails, which are paved with stones between sites at elevations of 
860-840 m (Fig. la). The trail impeded the movement of the soil to lower elevations from an elevation of 860 m, whereas the 
concentrated water flowing on the trail intensified the scouring of the soil at an elevation of 840 m. This is similar to a previous 
report [32] that the concentrated water flow of rural roads may scour the soil and affect the distribution of !?7Cs and the soil erosion 
rate. 
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Fig. 5 Distribution of slope gradient in the adjacent surface soil sample site in Longji Rice Terraces. (a) PZTF. (b) JRYTF. 


Fig. 5b shows the slope gradient distribution between two adjacent sites in the JRYTF. The slope first decreased and then 
increased with elevation from 905 m to 810 m, which is the opposite of the distribution trends shown for *°°*74°Pu and 137Cs 
dependence on elevation. The sampling site at an elevation of 1070 m (DA8) was close to the observation platform. Due to intense 
precipitation, irrigation, and other factors, surface runoff is accelerated by gravity and flows downward from the upper terraces, 
which greatly increases the runoff erosive power [56], resulting in lower activities of 7°°*°Pu and !37Cs at the sampling site (DA8). 
From elevations of 1070 m to 905 m (DA2), the activities of °9*?*°Pu and !37Cs increased with a decrease in elevation, which may 
have been caused by the gradual decrease in the erosivity of the surface runoff due to the gradual increase in the soil carried by the 
surface runoff. 


3.1.2. Vertical distributions of 7°°*?4°Pu and '°"Cs activities in soil cores 


The vertical distributions of °°*?4°Pu and '3’Cs activities in the PZTF soil core (P1 #1) at an elevation of 930 m and the soil core 
(P1 #2) at an elevation of 860 m are shown in Fig. 6a and 6b, respectively. The activities of 7°**4°Pu and '°’Cs in soil core P1#1 
tended to be uniformly distributed within the upper 20 cm and declined exponentially below this depth. The distribution of !37Cs 
activity in soil core P1#2 also tended to be uniformly distributed within the upper 20 cm, but the **°'*4°Pu value fluctuated greatly 
within the upper 20 cm, and both declined exponentially below this depth. The difference in the distribution of Pu isotopes and 
137Cs activities in the upper 20 cm layers of soil core P1#2 may be due to rainwater immersion. At site P1#2, water immersion 
enhanced the reducibility of the soil, leading to a change in Pu (IV) to Pu (IID, and the effect of oxide metals acting on Pu migration 
was weakened [57-59]. However, this soil environment had a very weak effect on '3’Cs because the reduction potential of Cs*/Cs 
was lower than that of the soil environment. Therefore, the difference in the distribution of Pu isotopes and '*’Cs in the upper 20 
cm layers of soil core P1#2 may be caused by the impact of rainwater immersion on the oxidation and reduction of Pu in the soil 
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environment. The oxidation and reduction of Pu in soil greatly influences Pu mobility [60]. The difference in the distribution of 
239+240Pu activity in the upper 20 cm layers from soil cores P1#1 and P1#2 may be due to discrepancies in the soil environment. 
P1#1 had fewer rainwater soaking times because of its location at the top of the terraces, which may have been less affected by 
rainwater immersion. Previous studies [21, 25, 55] have also shown that the vertical distributions of 73°*74°Pu and !37Cs in the plow 
layer were basically uniform. Therefore, the vertical distributions of 73°*74°Pu and !°7Cs in the two soil cores indicated that the 
terraced plow layer depth was approximately 20 cm, and 7°°***°Pu and !3’Cs were mixed uniformly within the upper 20 cm by tillage. 
The distribution of 73°*?4°Pu and !37Cs activities below the plow layer reflects their downward migration. The relatively high 7°*?*°Pu 
and !3’Cs activities within the upper 20 cm layer are also consistent with the '°’Cs results of the plow layer in the earth-banked 
terraces of the Yimeng Mountains [56]. 

The mean activities of °°*74°Pu and '3’Cs in the soil core P1#2 (0.242+0.252 and 3.58+2.17 mBq/g, respectively) are twice those 
of P1#1 (0.12140.080 and 1.84+1.46 mBq/g, respectively), and the maximum activities of °*°Pu and !3’Cs (0.69340.046 and 
5.77+0.24 mBq/g, respectively) are also higher than that of P1#1 (0.228+0.008 and 4.07+0.20 mBq/g, respectively), indicating a 
different extent of soil erosion at two sites. 


239+240py activity(mBq/g) 239+240py activity(mBq/g) 
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Fig. 6 Vertical distributions of %°'4°Pu and '°’Cs activities in soil cores from Longji Rice Terraces. (a) P1#1. (b) P1#2. 


3.2. Inventories of 2°°*24°Pu and !3’Cs in soil cores 


Calculated the inventories by integrating 7°°**4°Pu content and !7Cs content in all layers from the cores P1#1 and P1#2 (Fig. 6), 
we calculated that the inventories of 77°'*4°Pu in soil cores P1#1 and P1#2 were 81.51 Bq/m? (surface layer accounted for 14.87%) 
and 119.28 Bq/m? (surface layer accounted for 11.10%), respectively. And the '°’Cs inventories in soil cores P1#1 and P1#2 were 
1232.86 Bq/m? ( surface layer accounted for 16.42%) and 1827.61 Bq/m? ( surface layer accounted for 19.16%), respectively. 

Table 2 lists the locations, elevations, precipitations, and inventories of °*°Pu and '3’Cs of Longji Rice Terraces and the 
previous work. The mean inventories of 7°°*74°Pu and !7Cs in Longji Rice Terraces were 100.40 Bq/m? and 1530.24 Bq/m’, 
respectively. The mean inventory of 7°*?4°Pu was higher than the result (63 Bq/m’) of forest soil in Guiyang of Guizhou [42] but 


Table 2 Locations, elevations, precipitations, and inventories of °*?4°pu and !37Cs of Longji Rice Terraces (LJRT) and previous work. 


Sample Location Elevation Precipitation ?3%+?40Pu inventory !37Cs inventory reference 

site (m) (mm/y) (Bq/m?) Bq/m? 

Guiyang 106°40' 10"E, 1349 1130 63 - [42] 
26°39'23"N 

Kaiyuan 103°26'48.37"E, - 904 - 859.29 [61] 
23°29'20.15"N 

Libo 107°37'E-108°18'E, 860-1020 1470 - 723.29 [62] 
25°7'N- 25°9'N 

LJRT 110° 7'E, 860-930 1900 100.40 1530.24 This study 
25° 45'N 


The '?’Cs inventories were decay-corrected to 2021. 


within the range (36-410 Bq/m?) in the 20 °N - 30 °N band summarized by Cao et al. [63]. The mean inventory of '3’Cs was higher 
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than that previous studies in Kaiyuan [61] and Libo [62], but within the range (414.42-3140.88 Bq/m?, decay-corrected to 2021) in 
the 20 °N - 30 °N band summarized by Tagami et al. [64]. The global inventory of plutonium deposited at a specific location may 
be related to precipitation, elevation, latitude, and other factors [65]. As shown in Table 2, the latitudes of Longji Rice Terraces, 
Guiyang, Kaiyuan, and Libo are similar, and the elevation of Longji Rice Terraces is lower than that of Guiyang but similar to that 
of Libo. The precipitation of Longji Rice Terraces(1900 mm/y) was higher than that of Guiyang(1130 mm/y) and Libo(1470 mm/y), 
and the high plutonium inventory and '*’Cs inventory of Longji Rice Terraces may be due to the relatively high precipitation. 


3.3. Soil erosion rate 


The soil erosion rate of cultivated land can be estimated by a simplified mass balance model [4, 66, 67]: 

N-1963 
Im = I;(1 - =) (1) 
where Im is the measured inventory of 73°*74°Pu or !37Cs in the eroded soil core (Bq/m7), I; is the background value of 73°*°Pu or 
137Cs inventory (Bq/m?) in the study area, H is the depth of the plow layer (cm), A is the erosion thickness (cm/a), and N is the 
sampling year. Thus, the soil erosion rate EF (t/(ha:a)) is given by the expression [4]: 

E=100-h-B (2) 
where B is the soil bulk density (g/cm?). 

A background inventory (Z4) of °'**°Pu or '37Cs in the study area is crucial for estimating soil erosion. The traditional approach 
is to take an inventory of an undisturbed and flat reference site, which represents the background inventory of the study area. 
However, this approach is not applicable to the Longyji Rice terraces. If site P1#1 is taken as a reference site, we will obtain the 
unreasonable result that there is only soil deposition in the terraces, except for the top site P1#1 and the bottom site, because the 
lowest activity concentrations of 7°°*?4°Pu and '3’Cs on the surface (Fig. 2a) appeared at 930 m (the top site P1#1) and 790 m (the 
bottom site). It is also inappropriate to use site P1#2 as a reference because it is impossible to determine whether the soil flowing 
into and out of site P1 #2 is in equilibrium. In addition, the inventories of 7°°***°Pu and '°’Cs were higher than those of nearby areas 
(Table 2) because of the relatively high precipitation in the study area; therefore, neither of the sites listed in table 2 is an ideal 
reference site. It is also difficult to find an ideal undisturbed and flat reference site in the study area because of its unique erosional 
environment (special topography and geomorphological characteristics). Moreover, when the surface soil was eroded, the depth of 
the plow layer was approximately constant owing to the same tillage method; therefore, the deeper soil under the original plow 
layer was turned over and mixed evenly with the upper soil, which diluted the 73?*?4°Pu and !37Cs in the soil of the original plow 
layer, resulting in the redistribution of 7°**°Pu and !3’Cs in the surface soil from the terraces. Therefore, the distributions of 77°*4°Pu 
and !°’Cs in the surface soils from the studied area of Longji Rice terraces can reflect soil erosion, and the mean surface activities 
of 73°*240Py and !37Cs can be used as the background values of surface activity in the studied area. Therefore, we estimated the 
background values of 3°*74°Pu and !7Cs inventories from the mean activities of °'°Pu and '°’Cs, respectively, in surface soil. 

The 7°°**4°Pu in the surface layers (0-5 cm) of soil cores P1#1 and P1#2 accounted for 14.87% and 11.10%, respectively, of the 
total soil core inventory, with a mean value of 12.98%. The !°’Cs in the surface layers of soil cores P1#1 and P1#2 accounted for 
19.15% and 16.42% of the total soil core inventory, respectively, with a mean value of 17.79%. The inventories of 7°'*4°Pu and 
137Cs in the surface layer were 14.67 and 303.92 Bq/m/”, respectively, which were derived from the mean activity of 7°**°Pu (0.253 
mBq/g) and !3’Cs (5.24 mBq/g) as well as the bulk density (1.16 g/cm?) in the surface samples from PZTF. Then the background 
values of 113.05 Bq/m? for 73°*”4°Pu and 1708.38 Bq/m? for !°’Cs were estimated by the proportion (12.98% for °'**°Pu and 17.79% 
for '3’Cs, respectively) of inventory in the surface layer. The background values obtained using this method were slightly different 
from the mean inventories of °*°Pu and '°’Cs from cores P1 #1 and P1 #2 in the study area. For *°**°Pu, the calculated 
background value was 113.05 Bq/m? and the mean inventory from P1#1 and P1#2 cores was 100.4 Bq/m?°, and the difference 
between them was only 11.6%. For '°’Cs, the difference between the calculated background value of 1708.38 Bq/m? and the mean 
inventory of 1530.24 Bq/m? was only 12.5%. Therefore, the background values of the study area estimated using our method are 
reasonable and credible. 

Soil erosion rates were calculated using 7°°**4°Pu and !3’Cs methods (Table 3). The result showed that the erosion rates calculated 
by the 73°*74°Pu method were13.04 t/(ha-a) for P1#1 site and -2.15t/(ha-a) for the P1#2 site, which is consistent with the values of 
13.01 t/(ha-a) and -2.70 t/(ha-a) by '3’Cs method. The erosion rates were negative, indicating net soil deposition at site P1 #2. 

P1#1 was near the observation platform, and the slope gradient between it and the observation platform was steeper, resulting in 
severe erosion at P1#1 and no soil replenishment at that location. Sampling site P1#2 was located on the contour at an elevation of 
860 m (near PH11). The slope gradient was approximately 18.4°, which was gentler than the slope gradient between the upper 
terraces at higher elevations. Therefore, the eroded soil from the upper terraces was easily deposited in P1#2. Soil erosion rates at 
sampling site P1#1 were much higher than those at P1#2, again confirming that local landform features (e.g., slope gradient variation) 
are the main factors influencing soil erosion in the Longji Rice Terraces. 


Table 3 Soil erosion rates estimated by °*?4°Pu and !37Cs methods. 
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Sampling sites Erosion rates (t-ha!-a') reference 


239+240Py method 137Cs method 

P1#1 13.04 13.01 This study 
P1#2 -2.15 -2.70 This study 
Terraces in the Yimeng - 27 [56] 
Mountains 
Terraces in Zhongxian - 56.47° [26] 
County, Chongqing 
Dapotou depression, Yunnan - 6.32 [61] 
Cultivated land of Pingguo in - 19.4? [68] 
Guangxi 
Cultivated land of Huanjiang = - 0.22 [69] 
in Guangxi 

è mean value. 


The mean erosion rates (i.e., the arithmetic mean of the erosion rates at sampling points P1#1 and P1#2) of Longji Rice Terraces 
estimated by the 7°**4°Pu method (5.44 t/(ha-a)) and the !3’Cs method (5.16 t/(ha-a)) are basically consistent with that of the Dapotou 
depression in Yunnan Province (6.32 t/(ha-a)) [61], but higher than that of the cultivated land in Huanjiang County, Guangxi (0.22 
t/(ha-a)) [69]. The Dapotou depression area (23° 36'48.04"-23° 37'28.10" N) and the Longji Rice terraces are both hillside cultivated 


lands with similar latitudes, this is a possible reason why the mean erosion rates of both are similar. Higher precipitation (1900 
mm/y for Longji Rice terraces and 1390 mm/y for Huanjiang County) may lead to higher soil erosion in the Longji Rice terraces 
than in Huanjiang County. In addition, the reported cultivated land in Huanjiang County has been partially reclaimed, and some of 
the cultivated land was converted from cultivated land to grassland in 2004. The diverse land-use types and high vegetation cover 
may be other reasons for the lower soil erosion rate in Huanjiang County. This may be due to the higher vegetation cover of 
grasslands compared to cultivated land, and the interlaced plant roots in grasslands can inhibit soil erosion [36]. Moreover, the mean 
soil erosion rate of Longji Rice Terraces is much lower of the steep sloping terraces (56.47 t/(ha-a)) in Zhongxian County of 
Chongqing [26], the earth-banked terraces (27 t/(ha-a)) in the Yimeng Mountains [56] and the cultivated land (19.4 t/(ha-a )) in 
Pingguo City of Guangxi [68]. The soil of the steep-slope terraces in Zhongxian County is purple with loose parent rock that easily 
disintegrates and has faster soil loss. The Longji Rice Terraces are known as a hundred-stage terrace system, whereas the Yimeng 
Mountain terrace study area is only a seven-stage terrace system, which may have led to a lower erosion rate of the Longji Rice 
terraces. This may be because the terraces changed the slope gradient of the original hillslope, flattening the terraced plots at all 
levels. The more steps there were, the shorter the slope length of each terrace, and the terrace fields could slow or prevent rapid 
surface runoff due to the interception of runoff water and a shorter slope length, which significantly reduced runoff and soil erosion 
rates [70]. In addition, Longji Rice Terraces are cultivated once a year, while in Yimeng Mountain, they are cultivated twice a year 
(wheat-corn rotation), and lower tillage frequency may be another reason for the lower erosion rate in Longji Rice terraces. It has 
been reported that there is a sinkhole in the cultivated land in Pingguo City, and obvious underground erosion is a possible reason 
why the soil erosion rate in this area is higher than that of the Longji Rice Terraces. 

In summary, the results show that the differences in soil erosion rates in the above-mentioned areas are mainly related to 
precipitation, landform features, soil structure and properties, vegetation coverage, crop type, and tillage frequency. 


4. Conclusions 


The activities of 7°°*74°Pu and !37Cs in the surface soil were 0.089-0.734 mBq/g (n=33) and 1.80-7.88 mBq/g (n=33), respectively. 
All PuPu atom ratios were in good agreement with the global fallout value, indicating that plutonium in Longji Rice Terrace 
soils originates from the global fallout. The inventories of 7°°**°Pu and !°7Cs were 81.51-119.28 Bq/m? and 1232.86-1827.61 Bq/m’, 
respectively. The erosion rates calculated using the **°***°Pu method and '*’Cs methods were consistent. The distributions of 
239+240Py and !37Cs and the soil erosion rates at the different sample sites were mainly caused by local landform features, especially 
the slope gradient. The erosion rate of the P1#1 site calculated by the 7°°**°Pu method (13.04 t/(ha-a)) and the '°’Cs method (13.01 
t/(ha-a)) showed that there was severe soil erosion and necessary measurements should be taken to reduce it. Plutonium, as a 
potential tracer, can replace '3’Cs in the future and plays an important role in soil erosion research. 
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